Static vs. dynamical mean field theory of Mott antiferromagnets 
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Studying the antiferromagnetic phase of the Hubbard model by dynamical mean field theory, 
we observe striking differences with static (Hartree-Fock) mean field: The Slater band is strongly 
renormalized and spectral weight is transferred to spin-polaron side bands. Already for intermediate 
values of the interaction U the overall bandwidth is larger than in Hartree-Fock, and the gap 
is considerably smaller. Such differences survive any renormalization of U . Our photoemission 
experiments for Cr-doped V2O3 show spectra qualitatively well described by dynamical mean field 
theory. 

PACS numbers: 71.27.+a, 71.20.-b, 79.60.-i 
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I. INTRODUCTION 

In recent years, our ability to calculate strongly corre- 
lated materials has substantially improved. To this end, 
one needs to go beyond the conventional local density ap- 
proximation (LDA) ii New methods had to be developed 
like 1jDA-\-U ^ where LDA is supplemented by a local 
Coulomb interaction U treated in the static Hartree-Fock 
(HF) mean field theory, and its sibling LDA-FDMFT^ 
which employs the more sophisticated dynamical mean 
field theory (DMFT)*^ While it is generally accepted that 
LDA+DMFT deals more accurately with strongly corre- 
lated metals, the simpler LDA-I-U is considered to be suf- 
ficient for insulators with a large U ^ at least in the pres- 
ence of magnetic or orbital ordering. Indeed. LDA-fU 
was tailored for such strongly correlated insulators^ and 
is nowadays widely employed to calculate various physi- 
cal quantities of these. 

In our paper, we work out the differences be- 
tween a static and a dynamical mean field treat- 
ment of long-range-ordered insulators at intermediate 
to strong Coulomb interactions U . For realistic mate- 
rial calculations, this would correspond to LDA-I-U and 
LDA+DMFT. We study these differences, considering 
the antiferromagnetic (AF) phase of a simple model, the 
one-band Hubbard model. We point out that at large 
U dynamic properties and the structure of excited states 
are strikingly different in the Slater (HF) and the DMFT 
antifcrromagnet. More specifically the Slater bands are 
strongly renormalized, most of the spectral weight is 
transferred to spin-polaron side bands, and the overall 
bandwidth is proportional to the non-interacting width 
W, as opposed to the 1/U shrinking found in HF. We 



also performed photoemission spectroscopy (PES) exper- 
iments for V2O3 doped with 1.2% Cr, considering both 
the AF and the paramagnetic insulator. These experi- 
ments confirm that strongly correlated antiferromagnets 
are not of Slater (HF) type, while they can in many re- 
spects be described by DMFT. 

The paper is organized as follows: In Sec. ^the DMFT 
results for the antiferromagnetic phase of the Hubbard 
model are presented and compared to the exact solution 
of the t-J model in infinite dimensions. In Sec. Illll we an- 
alyze the evolution of staggered magnetization and spec- 
tral function from weak to strong coupling. This allows 
us to draw some general conclusions about the validity of 
DMFT and HF in describing both ground- and excited- 
state properties in the antiferromagnetic phase. In Sec. 

the evolution from a Mott antifcrromagnet to a para- 
magnetic Mott insulator with increasing temperature is 
discussed. The photoemission spectrum of V2O3 is also 
presented in both phases and it is compared to our theo- 
retical calculation. Finally we present our conclusions in 
Sec. M 



II. ANTIFERROMAGNETIC PHASE OF THE 
HUBBARD MODEL 



Several studies have been performed on the AF phase 
of the Hubbard model|^ but here we focus on large values 
of U , that received less attention in the past, and we 
exploit some technical advances to improve the accuracy 
of our DMFT calculation so that we can make definite 
statements not only about the size of the AF gap but also 
regarding the inner structure of the (Hubbard or Slater) 
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exact t-J : J=D/8 

2 - DMFT ; U=8D ; N =25 
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FIG. 1: (Color online) Spectral function of the lower Hubbard 
band. Red (dark gray) curve: DMFT solution of Hubbard 
model for U = 87? and Ns = 25. Green (light gray): t-J 
model for J = /U. Both spectra have been plotted using 
the same Lorentzian broadening. 



bands below (and above) the gap. In the large U limit, 
we can make contact with the t-J model which has been 
studied in infinite dimensions by Strack and VoUhardti 
We recover the spin-polaron peaks of Ref. which are 
however dispersive in our calculation, an effect occurring 
in order D'^/U^ and hence absent in Ref. 0. 
The Hubbard Hamiltonian reads 

H^-t 4aCja + H.C. + UY, (1) 

where t is the hopping amplitude, U is the Coulomb re- 
pulsion, Cia (cjg.) are annihilation (creation) operators 
for spin cr electrons on site i and rii^ = cl^Cia- We solve 
the DMFT equations for a semicircular density of states 
N{e) = {2/TTD^) V£)2 _ £2 ^iti^ bandwidth W = 2D, us- 
ing exact diagonalization and Lanczos algorithm at T = 
for the associated Anderson impurity model (AIM). A 
continued fraction expansion of the DMFT Weiss fields^ 
allows us to obtain reliable spectra in the intermediate- 
to-large U region. Moreover, in the AF phase we are able 
push the number of AIM sites Ng to much larger values 
than for the paramagnetic phase. 

To confirm that we get reliable spectral functions we 
first consider the laige-U limit of the Hubbard model 
and compare it to the exact solution for the t-J model in 
infinite dimensionsii The comparison of the fc-integrated 
spectral functions is displayed in Fig.^ In both cases, we 
show only the negative frequency part of the spectrum, 
i.e., the lower Hubbard band centered around uj = —U/2 
(the X axis is shifted correspondingly) . Evidently in both 
cases we find spectra of the same type. We have also 
verified that the DMFT spectra are almost independent 
on Ns, reflecting the fact that, contrary to the case of 
the paramagnetic insulating phase, in the AF phase the 



peaks do not originate from the discreteness of the im- 
purity model but have a physical meaningiS 

As already observed in Ref. they originate from the 
fact that a hole moving in the ordered background breaks 
AF bonds costing an energy proportional to nJ, where n 
is the length of path in the lattice. This string potential 
gives rise to a set of discrete energy levels with a typical 
separation proportional to J^/'^. These levels can be in- 
terpreted as spin-polaron side peaksifi for dispersionless 
spin-waves^ In the U oo (J 0) limit, these spin- 
polaron peaks become dense and the lower Hubbard band 
recovers the shape and the width of the non-interacting 
density of states, with the important difference that the 
states are incoherent for [/ — > cxd. In our faiite-U case, 
we find that the shift of the peaks in the Hubbard model 
with respect to the t-J model scales with l/C/^. That 
means we recover the infinite-?/ limit where the mapping 
of the Hubbard to the t-J model is exact. 



III. STATIC AND DYNAMICAL PROPERTIES 
FROM WEAK TO STRONG-COUPLING 

After having demonstrated the accuracy of our spectra 
we investigate how the picture changes when U is lower 
and follow the evolution of our antiferromagnetic solution 
as a function of U. 

While, for f7 = D static and dynamical mean-field the- 
ories yield, as one expects, similar spectra (see uppermost 
curve in the main panel of Fig. [SJ, already for U = 3D 
the spin-polaron picture characteristic of the t-J model is 
almost fully developed. This is instead completely missed 
by the static mean-field theory. The discrepancy between 
DMFT and HF spectrum grows with increasing U, as ex- 
emplified in the figure for U = 8D. The HF spectrum 
arises from a renormalized single-particle band with a 
dispersion decreasing as 1/U. In contrast in DMFT the 
width of the Hubbard band is always of the order of the 
bare bandwidth even in the limit U ^ oo. 

Focusing on the rightmost peak of the Hubbard band, 
we find that it continuously evolves from the weak- 
coupling Slater (HF) peak. However, its spectral weight 
Z, directly measured from the spectral density, is dra- 
matically reduced when increasing U since more and 
more weight is transferred to the spin-polaron side bands. 
As it can be seen in the inset of Fig. |21 Z scales as D/U 
for U > 2D. This is what is also found for the t-J model, 
however the physical nature of the excitation is different. 
In the t-J model the peaks are dispersionless^ while in 
the Hubbard model we find, for the peak closest to the 
Fermi level, a bandwidth w which scales as D^/U"^ for 
large U (solid diamonds in the inset of Fig. ^ . 

An independent argument for the dispersion of this 
peak is obtained by explicitly looking for the poles of the 
Green's functioni^ Averaging over the two sublattices we 
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FIG. 2: (Color online) Evolution of the lower Hubbard band 
with U/D. The solid gray lines represent the HF spectrum 
for U = D and U = 8D. In the inset the spectral weight and 
the quasiparticle dispersion are shown as a function of U/D 




FIG. 3: (Color online) DMFT spectral gap as a function of 
the AF order parameter Ao = Um, where m is the staggered 
magnetization. In the Hartree-Fock case the gap coincides 
with Ao (dashed line). In the inset we compare the DMFT 
and the HF magnetization as a function of U. 
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where is the bare band energy, = cj + /.t — Sq(w) 
and a = A, i? is a sublattice index. The renormalized 
dispersion E-^ is given by the solutions of 



Rea(^)CB('^) 







(3) 



with Ek ranging from ~D to D. 

We label the solution for = 0. At large U we 
expand Eq. 10) to linear order and obtain 



Ek 



where we have defined 
d 



du 



ReCA(w)Ci3('^) 



(4) 



(5) 



To determine the value of >c, we consider the Green's 
function for energies ui ~ —U/2, which for large U is 
given by the Green's function of a given sublattice. e.g. 
Gaa{uj). In this limit, Eb « -Sa ~ U/2. Then we find 
that the quasiparticle weight is given by Z = —U/k. 

In particular, considering as an example the case of a 
hypercubic lattice, the dispersion E^ has the following 
form: 



E-k — const. 



(6) 



where fc„ is the wave vector in the nth dimension. 



The presence of a finite quasiparticle dispersion in 
the Hubbard model can be understood as an effect of 
the 0(1/?/^) terms neglected in the mapping to the t-J 
Hamiltonianpi^ More precisely, in addition to the stan- 
dard t-J terms, the large-C/ projection of the Hubbard 
model leads to spin-flip terms and to second (^2) and 
third (t-i) nearest neighbor hopping integrals. Such high- 
order hoppings contribute to the quasiparticle dispersion 
even in infinite dimensions, and in the specific case of a 
hole in a Neel state read = 2t^ /U and ^3 = t'^ /U, re- 
spectively. For a hypercubic lattice these terms give rise 
to a bare dispersion of the form 



-2t2 cos kn COS km — 2t^ cos 2k„ 



n^ra 



= — — (2f cosfc„)^ — const. 



(7) 



From the results obtained in Rcf. for a t-t'-J model, 
we expect for the Hubbard model, in the limit of t2 <C 2J 
and ts ^ 2J, a dispersion given by the bare of Eq. 
{Tj) renormalized by the quasiparticle weight Z ~ D/U . 
This expectation is in fact confirmed since Zik coincides 
with the dispersion given by Eq. JB)). 

These considerations support the interpretation of our 
DMFT results as a quasiparticle renormalization of the 
Slater (HF) peak; the width of this peak shrinks like 

/U and it is further reduced by Z ~ D/U so that 
altogether w~ZD'^/U = D^/U'^. 

As already mentioned before, a striking difference be- 
tween the Slater and the DMFT spectrum is the ampli- 
tude of the spectral gap. Since the HF Slater peak is 
located roughly at the center of the DMFT lower Hub- 
bard band of width of order W , for large U the spectral 
gap is reduced from the HF value U toU — W . We illus- 
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trate this aspect in Fig. where the amplitude of the 
spectral gap is plotted as a function of the order param- 
eter Aq. Aq is given by U multiplied by the staggered 
magnetization to, whose behavior is shown in the inset 
of Fig. O The DMFT results indicate that, on the one 
hand the staggered magnetization gets closer and closer 
to the Hartrec-Fock value by increasing U, while, on the 
other hand, the spectral gap deviates more and more 
from the Hartree-Fock prediction (Aq) approaching the 
asymptotic value of Aq — W. Thus we have the interest- 
ing situation where the two methods give essentially the 
same magnetization, while there is a sizeable difference 
in the gaps. 

Our findings highlight one important point: At large 
U, DMFT and HF give the same description of the anti- 
ferromagnetic ground-state, which is basically the Neel 
state. On the other hand, the description of excited 
states is strikingly different between DMFT and HF. In 
HF the hole moves in a "rigid" background of Neel spins, 
while, in DMFT, due to quantum fluctuations, the exci- 
tations are of a completely different nature: The hole can 
also move around without paying the double occupancy 
cost, with a creation of a spin polaron. Therefore, even if 
all ground-state quantities (to, kinetic and potential en- 
ergy, . . . ) are basically the same in the two descriptions, 
excited state properties (hence spectra) are completely 
misrepresented by HF. 

This also means that, within HF, it is not possible 
to simultaneously obtain the correct magnetization and 
spectral gap with the same value of U. This has obvious 
implications for the use of LDA-I-U for antiferromagnets 
(or similarly for orbitally-ordered phases) with a large 
gap. 



IV. ANTIFERROMAGNETIC VS 
PARAMAGNETIC MOTT INSULATOR 

Our DMFT calculations can be extended also to fi- 
nite temperature. In this case, calculations can only be 
performed with a lower value of iV,, but their precision 
is still sufficient to address a crucial issue of the Mott 
antiferromagnet, namely the evolution from a, t-J anti- 
ferromagnet to a paramagnetic Mott insulator above the 
Neel temperature. In particular we are interested in the 
mechanism ruling the evolution of the system from a fully 
polarized antiferromagnet with a spectral gap of order U 
to a Mott insulator with zero magnetization but with a 
gap of the same order of magnitude. At low T, due to 
the onset of long-range AF order, the opening of the gap 
is accompanied by a rigid shift in the real parts of the 
self-energy, whose value is given by the magnitude of the 
order parameter. The imaginary parts are, in this case, 
small for both spin species and non-zero only where the 
corresponding (up or down) spectral weight is present 
(upper panel of Fig.^J. On the other hand, for T > Tjy, 
when the order parameter is zero, a Mott gap of order 
U is found, associated with a ^ U'^ /lu peak in the imag- 
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FIG. 4: (Color online) Imaginary part of up and down self- 
energy for U = 5D, below (two upper panels) and above (low- 
est panel) the Neel temperature {Tm = 0.05D). The spectral 
gap is marked by the dashed vertical lines and Aq — Um. 



inary parts of T,{uj) at zero frequency, as can be seen in 
the lowest panel of Fig. 21 How can we go from one sit- 
uation to the other? By increasing T, two peaks in the 
imaginary parts of S((x') develop inside the spectral gap 
(indicated by the vertical dashed lines in the figure) and 
their mutual distance decreases together with the order 
parameter Aq: The two peaks can be viewed as precur- 
sor of the ~ U'^/u! behavior in the paramagnetic Mott 
insulator. 

Let us now ask: Are real transition metal oxides de- 
scribed by static or dynamical mean field theory? To 
this end, we performed photoemission measurements for 
V2O3 doped with 1.2% Crii4 The measurement was 
made at beamline BL25SU of SPring-8 using 500eV pho- 
tons. The experimental details are exactly as described 
previouslyii^ The comparison between theory and exper- 
iments is shown in Fig. taking parameters reasonable 
for V2O3: U = 3[eV] and 2D = 2[eV]. As we discuss be- 
low, the theoretical spectra have been broadened to make 
contact with the experiment. Both theoretical and exper- 
imental spectra are normalized to have the same area. A 
Shirley-type inelastic background has been removed from 
the experimental data. 

Comparing the experimental data (right panel of 
Fig.©) to the spectra we have discussed until now, it is 
evident that the spiky nature of the latter is absent, or at 
least basically invisible, in the former. This is evidently 
due to at least three sources of broadening: the intrinsic 
experimental broadening, the finite temperature and fi- 
nite dimensionality effects beyond the DMFT. This last 
source of broadening can, e.g., arise from the coupling to 
dispersive spin- waves ^ The size of this effect is hard to 
estimate for three-dimensional compounds. However, to 
mimic the overall effect we have plotted in the left panel 
of Fig.© our theoretical results using a broadening of 
0.2[cV], larger than the experimental one (~ 0.1[eV]) in 
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FIG. 5: (Color online) Comparison of theory (left) and PES 
(right) spectra of the lower Hubbard band in AF and para- 
magnetic phases. The black, cyan (light gray), and magenta 
(dark gray) arrows, indicate, respectively, the formation of 
the (quasiparticle-)renormalized Slater peak, the reduction of 
the gap in the AF phase and the qualitative disagreement in 
the lower edge of the Hubbard band. 



order to include phenomenologically the other effects. 

A first observation is that the experimental spectrum 
in the antiferromagnetic phase has an overall width of the 
same order of the paramagnet. We have already shown 
how this is realized in the DMFT spectra, in contrast to 
the static HP mean-field. Theory and experiment qual- 
itatively agree also on the substantial redistribution of 
spectral weight: Going from the paramagnetic to the AF 
insulator, an additional peak develops at the upper edge 
of the lower Hubbard band. This peak is the one we 
identified as a (quasipartiele-) renormalized Slater band. 
Besides this peak there is additional spectral weight com- 
ing from the broadened spin polaron side-bands. Re- 
markably, the gap is found to increase when going from 
the paramagnetic to the AF insulator both in experi- 
mental and theoretical spectra, ft is worth remarking 
that this feature does not depend on a particular choice 
of the value of U , since it occurs even in the infinite-t/ 
caseii^ The agreement is less good in one qualitative as- 
pect: the flank at the lower edge of the Hubbard band is 
almost identical for both phases in experiment whereas 
they slightly differ in the theoretical calculation. This is 
likely due to the lack of some realistic features, like an 
asymmetric density of states and multi-orbital effects in 
our model calculations. This point could be clarified by 
extending LDA-I-DMFT calculations from the paramag- 
netic phasoii to AF V2O3. 



V. CONCLUSIONS 

In this work we have compared static and dynamical 
mean field descriptions of the antiferromagnetic phase 



of the Hubbard model. We have analyzed the evolution 
from weak to strong coupling of both ground-state (e.g. 
magnetization) and excited-states (e.g. spectra) proper- 
ties. Contrary to the general expectation that HF gives 
a good description of broken-symmetry phases both at 
weak and strong coupling, we have found that, at large U, 
this is true only for ground-state quantities. The poor de- 
scription of the excitations given by static-mean field re- 
flects in a HF spectrum given only by two unrenormalized 
{Z = 1) Slater bands which shrink as 1/U with increas- 
ing U. On the contrary, when the excited states are bet- 
ter described, as in DMFT, the Slater bands get strongly 
renormalized {Z ~ D/U) and most of the spectral weight 
is transferred to spin polaron side bands. Therefore the 
total width of the Hubbard band stays finite (of order 
W) even for large U and the spectral gap is not given by 
U , as in HF, rather by [/ — W. 

We have also performed PES experiments on Cr-doped 
V2O3 both in the antiferromagnetic and paramagnetic 
phase. We have found a qualitative agreement with the 
DMFT calculations, in the formation of the renormalized 
Slater peak, in the overall width of the band (of order W 
rather than D/U), and also in the shrinking of the gap 
when going from the paramagnetic to the antiferromag- 
netic phase. 

In conclusion, we demonstrated with our DMFT cal- 
culation and by comparing theory with experiments that 
one has to be careful when applying static mean field the- 
ories, like LDA+[/, for calculating physical properties of 
insulators with long range order. Adjusting U to get the 
correct gap, one obtains wrong estimates of the magneti- 
zation, or the energy, and vice versa. These problems can 
be overcome by employing dynamical mean field theory. 
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